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Abstract 

We study finite temperature effects in two distinct holographic models that 
exhibit Lifshitz scaling, looking to identify model independent features in the 
dual strong coupling physics. We consider the thermodynamics of black branes 
and find different low-temperature behavior of the specific heat. Deformation 
away from criticality leads to non-trivial temperature dependence of correlation 
functions and we study how the characteristic length scale in the two point 
function of scalar operators varies as a function of temperature and deformation 
parameters. 



1 Introduction 

A number of holographic models have been introduced to study quantum critical sys- 
tems in low dimensions, see [IHH] for reviews. In the holographic approach, relatively 
simple classical or semi-classical calculations are carried out in a suitably chosen grav- 
itational background and the results interpreted in terms of strongly coupled physics 
in the dual field theory. For condensed matter applications, the holographic models 
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are usually of the bottom up variety, where key symmetries and conserved charges of 
the dual system under study are realized by coupling gravity to a minimal number of 
additional fields. Without the benefit of extended supersymmetry and the pedigree 
of a consistent embedding into string theory, the validity of the gauge theory /gravity 
correspondence is very much an open question for bottom up models. On the other 
hand, they offer a simple phenomenological setting in which to study a range of phys- 
ical effects and one hopes to identify a set of generic features in bottom up models 
that may also be shared by more realistic top down constructions. Much of the work 
in this area has focused on systems with an underlying conformal symmetry. This 
is partly because conformal symmetry is realized in a number of interesting physi- 
cal systems, but also because the assumption of conformal symmetry often simplifies 
calculations. On the other hand, conformal symmetry is usually not present in con- 
densed matter systems and generic quantum critical points exhibit scale invariance 
without conformal invariance. 

In a quantum phase transition, the ground state of a system changes abruptly 
when some parameter of the system is varied through a critical value (for reviews, see 
e.g. [HIE])- As in a conventional phase transition, characteristic length scales of the 
system diverge at a quantum critical point and there is an emergent scaling symmetry 
of the form 



where z > 1 is referred to as the dynamical critical exponent. For z > 1 the scaling at 
the quantum critical point is anisotropic between the temporal and spatial directions, 
commonly referred to as Lifshitz scaling in the literature. 

In this paper we study finite temperature effects in two different bottom up holo- 
graphic models that exhibit Lifshitz scaling. On the one hand we consider Einstein- 
Maxwell gravity coupled to a massive vector (Proca) field. We refer to this as the 
EMP model (for Einstein-Maxwell-Proca). It is based on a model introduced in [8] 
involving a pair of two- and three- form field strengths coupled to Einstein gravity with 
a negative cosmological constant. This was reformulated in [9] as Einstein gravity 
coupled to a Proca field, and in [TU] a U(l) Maxwell field was added in order to study 
dual systems at finite charge density. The second model was introduced in [TT] and 
consists of Einstein gravity coupled to a dilaton field and a pair of U(l) gauge fields. 
It is a generalization of an earlier model considered in [9] and will be referred to below 
as the EDM model (for Einstein-Dilaton-Maxwell). A motivation for studying two 
models at once is to get a handle on how much of the physics is controlled by the 
anisotropic scaling invariance ([T]) and how much is model dependent. 

We compare charged black brane solutions and some of their thermodynamics be- 
tween the two models. We then consider two point correlation functions of operators 
dual to a probe scalar field in these holographic models. The two point correlators 
exhibit screening at finite temperature and we obtain the dependence of the screening 
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length on temperature and deformation parameters in the two holographic models. 
These results can be compared at a qualitative level to corresponding computations 
at weak coupling. In particular, the quantum Lifshitz model provides a well con- 
trolled theoretical setting to study thermal effects near a. z = 2 fixed point [12HT3]. 
The temperature dependence of characteristic length scales was studied via the two 
point function of a scaling operator in this model in [15J. It was found that the two 
point correlation function at spatially separated points vanishes at finite temperature 
due to an infrared divergence. This makes the system very sensitive to deformations 
away from the critical point. The authors of [15] considered the effect on the two 
point function of introducing an irrelevant operator and found that it restores corre- 
lations between spacelike separated points and introduces a temperature dependent 
characteristic length scale into the two point function. 

In general, heating up a system at a quantum critical point introduces a length 
scale that breaks the scaling symmetry, but at the same time it broadens the region 
in parameter space that is controlled by the underlying scale invariant theory. In an 
otherwise scaling symmetric theory, the dependence of a characteristic length scale £ 
on temperature is given by 

£ = cT~ l '\ (2) 

where c is some constant numerical coefficient. As one deforms the system away 
from the critical point, the temperature dependence of the characteristic length scale 
becomes more generic. Consider a set of deformation parameters Aj, which for conve- 
nience are defined so that they have dimensions of inverse length. At zero temperature 
Aj = corresponds to infinite correlation length. Away from A« = 0, the temperature 
dependence becomes 

Z = T- x '*7i{T- x t*\), (3) 

where rj is some generic function of its arguments. The only requirement is that 
77(0) = c. 

In the holographic models, turning on a finite temperature leads to the exponential 
decay of two point correlation functions between spacelike separated points, with a 
temperature dependent correlation length £. If the finite temperature is the only 
deformation parameter, then the temperature dependence of the correlation length 
will be of the simple form (J2}. More interesting behavior is seen if we consider 
deformations away from the fixed point. We study two different types of relevant 
deformations of the model. The first one is to turn on a chemical potential for a 
conserved U(l) charge, 

Si = J dxdrup, (4) 

which on the gravity side corresponds to considering gravitational solutions with an 
electric flux. The relevant black brane solutions in the EMP model were constructed 
in [10] and in the EDM model in [TTJ . As the temperature is lowered, this deformation 
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will take the system from the Lifshitz fixed point to a "locally critical" fixed point 
corresponding to the appearance of an AdS% x R 2 region in the spacetime. 

The second type of deformation we consider is to add a double trace deformation 
for the probe scalar field 

S 2 = J dxdrXO 2 , (5) 

whose correlation functions we are interested in. This deformation is less dramatic 
than the previous one in that, as the temperature is lowered, the theory flows to 
another Lifshitz symmetric fixed point. 

The paper is organized as follows. In section [2] we briefly review the calculation 
of the finite temperature correlation function in the quantum Lifshitz model in |15j . 
In section [3] we introduce the two holographic models with Lifshitz scaling and re- 
view properties of their charged black brane solutions. In section H] we compare some 
thermodynamic properties of black branes and conclude that the low temperature 
physics differs between the two models. In section [5J we obtain two point correlation 
functions in a charged black brane background. We first consider operators of large 
scaling dimensions, for which the two point correlator can be obtained in a geodesic 
approximation, and then repeat the analysis for operators of general scaling dimen- 
sions. In section Owe calculate the two point correlation function with a double trace 
deformation at finite temperature (and at a vanishing charge density). Finally in 
section [7| we summarize our results. 



2 Thermal correlators in the quantum Lifshitz model 

A free field theory realization of a class of z = 2 quantum critical points is provided 
by the quantum Lifshitz model, 

S = \ J d 2 xdr{{d T x) 2 + K(V 2 X f). (6) 

Here we will briefly review the finite temperature correlation functions in this model, 
obtained in [15], in order to compare them to the results of holographic calculations 
later in the paper. We refer to [T2TTT4] for more detailed discussions of the quantum 
Lifshitz model. 

The operators of interest are the so called monopole operators^ 

0(x) = e 2 ™ xix) . (7) 
The vacuum two point correlation function has the scaling form 



G(x,x') = (C(x,0)C(x / ,0)) oc 




1 Thc name monopole operator comes from a dual gauge theory representation of the same model, 
where one has Ei — tijdjX- The operator in (f7|) creates a monopole in the dual gauge field [15] . 
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which follows from logarithmic correlations of X- At finite temperature the contribu- 
tion from the zeroth Matsubara mode to the x correlator is proportional to 



T / d 2 k 



I _ gik-(x-x') 



(9) 



where the momentum integral is seen to diverge logarithmically with an infrared cutoff 
L. The two point correlation function of spacelike separated monopole operators thus 
vanishes at finite temperature in the infinite volume limit, 



G(x,x') oc 



e -f£l x - x 'l 21 



-> 0. 



(10) 



Spatial correlations are restored away from the scaling limit of the free theory. The 
authors of [15] consider a marginally irrelevant deformation of the free model, by 
adding an interaction term 



which corresponds to slightly moving away from the critical point. This deforma- 
tion leads to a logarithmic violation of scaling symmetry, which changes the finite 
temperature behavior. The interaction induces a term 



v\V X f 



(12) 



into the effective action for x already at one loop. This term changes the low mo- 
mentum behavior of the theory dramatically and cures the infrared divergences. By 
using a large N gap equation the authors of pE] find 



v 2 



log - logT 



(13) 



-logT 

which is understood to hold for small T. This means that there is a non-trivial 
(temperature dependent) length scale induced in the theory 



e-- 



(14) 



Using the new effective action leads to a two point correlator of monopole operators 
of the following form: 



G(x, x') oc exp 



7T 



T |x 



2v 2 



log 



2fr 



+ c 



where c is a constant, for distances smaller than £y, and 



G(x, x') oc exp 



2ttT 

— 5~ lo g 

v 



(15) 



(16) 
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for distances larger than Furthermore, if one allows for vortices (spinons) in the x 
field, the correlation function changes qualitatively at very large distances. Assuming 
a vortex energy gap E c , one finds that correlations are exponentially decaying on a 
length scale £, vor tex oc e Ec ^ 2T due to a vortex plasma [15] . 

To summarize, the authors of [15J find a non-trivial temperature dependence in 
correlation functions of spatially separated monopole operators, which arises on the 
one hand from the Lifshitz scaling region (at short distances) and on the other hand 
from deformations away from criticality. It seems likely that many of the details 
depend on the specific theory and the specific deformations that were considered 
in [15]. In sections |5] and [6] below, we calculate two point correlation functions of 
scaling operators in holographic models exhibiting Lifshitz scaling. We also find a 
temperature dependent thermal length scale which can be compared to that of the 
quantum Lifshitz model. Since the holographic models are dual to strongly coupled 
theories and the quantum Lifshitz model is a free field theory, one does not expect to 
reproduce the results of [15] in detail but rather look whether the holographic models 
exhibit similar trends. 



3 Holographic models with Lifshitz scaling 

We will consider two different holographic models, both of which realize anisotropic 
scaling of the form ([1]) through spacetime geometries that are asymptotic to the so- 
called Lifshitz geometry [5|ITB]. 

ds 2 = e (-r 2z dt 2 + ^ + r 2 dx 2 ) . (17) 

Here £ is a characteristic length scale of the geometry and the r, t, and x coordinates 
have no length dimensions. The Lifshitz metric is invariant under the transformation 

T 

t^\ z t, x — >■ Xx, r ->• —, (18) 

A 

which incorporates the scaling in ([T]) on the t and x coordinates, while the radial 
coordinate r of the bulk geometry scales inversely with A. In the rest of the paper 
the characteristic length will be set to I = 1. It can be re- introduced via dimensional 
analysis if needed. For concreteness, we will take the bulk spacetime to be 3+1 di- 
mensional, but our results generalize in a straightforward fashion to other dimensions. 

The two models both involve Einstein-Maxwell gravity with a negative cosmolog- 
ical constant but the asymptotic Lifshitz behavior is achieved by coupling to different 
matter sectors. Both models have charged black brane solutions that are asymptotic 
to the Lifshitz geometry and are interpreted as holographic duals of field theory con- 
figurations at finite temperature and finite charge density. The black brane geometries 
are not the same in the two models. In particular, their near-extremal limits differ and 
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this has consequences for low-temperature physics in the corresponding dual field the- 
ories. One of the aims of this paper is to compare predictions for physical observables 
at finite temperature and explore to what extent they are model dependent. 

3.1 Einstein-Maxwell-Proca theory 

The first model we consider is a modified version of the holographic model first con- 
sidered in [8] . In the original version, anisotropic scaling of the form (ITS]) was obtained 
by coupling Einstein gravity to a pair of two- and three-form field strengths. As was 
shown in [9], this can be rewritten as a single massive vector field coupled to gravity 
and we find this formulation more convenient to work with. Finally, following [10J, 
we include a Maxwell gauge field which is coupled to the gravitational field but not 
directly to the auxiliary massive vector field, 

Semp = j d 4 x^ (r - 2A - \f, u F^ - \j^ v - ■ (19) 

The role of the massive vector field background is to modify the asymptotic behavior 
of the metric. One could in principle include direct couplings between the massive 
vector and the Maxwell gauge field, in which case having a massive vector background 
field would also influence the gauge field dynamics, but for simplicity we have chosen 
not to do that. In what follows we will refer to this theory as the EMP model. 

The Lifshitz metric (ITTj) is a solution of the equations of motion when the constants 
in the model are related to the dynamical critical exponent z through 

A = - X -{z 2 + z + 4), c=V2i, (20) 
and the massive vector field has a non-vanishing background value, 

A t = ^^r z , A Xl =A r = 0. (21) 

The field equations of the EMP model have been studied extensively by many authors 
(see [17] for a brief review) and we will not repeat the analysis here. For generic values 
of the dynamical critical exponent, there are no known analytic solutions describing 
black holes or black branesj but numerical solutions can easily be found using similar 
techniques as introduced in [19]. Below, we will present results obtained from nu- 
merical solutions of the EMP model and compare them with results obtained using 
analytic black brane solutions from the other model. 

2 An exact solution exists for z = 4 (or more generally z = 2d where d is the number of transverse 
spatial dimensions) and a particular value of the electric charge on the black hole, or charge density 
on the black brane [lu]H8]. 
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3.2 Einstein-Dilaton-Maxwell theory 

The second model is an Einstein-Dilaton-Maxwell (EDM) theory 

2 

I 



EDM 



d 4 Xy/—g 



R-2A- -d u 
2 M 



1 2 



1=1 



(221 



Models of this type were introduced in the context of non-relativistic holography 
in [9] and black brane solutions of this particular action were found in Dynam- 
ical solutions describing infalling energy were constructed in [20] and used to study 
thermalization following a non-relativistic holographic quench but here we will only 
consider static solutions. 

The equations of motion are 



A,; 



A 9, 



V 2 - T e ^ F tt Fi 



1=1 



_i_ t^C 1 ) _i_ 7^( 2 ) 
' ~ //f ' 



0. 

0, 2 = 1,2, 



with 



T, 



fJ,V 



1 1 



rp(i) 



_ n p(i)p(i)crp 



(23) 

(24) 
(25) 

(26) 
(27) 



Remarkably, these field equations can be solved analytically to obtain charged black 
branes in asymptotically Lifshitz spacetime for any value of the dynamical critical 
exponent z. Below we reproduce the static black brane solutions of [11] and review 
some of their thermodynamic properties. We then use them to calculate finite tem- 
perature correlation functions of scalar operators in the dual field theory. At each 
stage, we compare our results to corresponding calculations in the EMP model. 



3.2.1 Asymptotically Lifshitz solutions 

We use the following metric ansatz for static solutions 

dr 2 



ds 2 = -r 2z f(r)dt 2 + 



r 2 g(r) 



+ r 2 (dx 2 + dy 2 



for which the field equations (1231) for the gauge fields are solved by 



r 



9{r) 



(28) 



(29) 
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where r is an arbitrary reference value of the radial variable and the pi are dimen- 
sionless constants. Under the coordinate rescaling fllSp . r transforms like r while pi 
remain invariant. Later on, when we consider black hole solutions it will be convenient 
to take ro as the radial location of the event horizon. 

The Einstein equations (123]) reduce to a pair of first order differential equations, 

y-y = -2(,-l) + ^(r0') 2 , (30) 



r f r Q' n/ ~\ 1 

J r + — = -2(z + 2)-- 

f 9 9 
and the dilaton equation (1241) becomes 



2K + \(^)\ple-^ + ple-^) 



(31) 



r 2 <b" + 



rf r g' 



+ Y g {j}' M" - *'* + = 0. (32) 

The Lifshitz geometry (fT7|) has f(r) = g(r) = 1, for which equation (|3"U1) is solved by 
e <t> = ^ r 2Vz-i_ fj^g integration constant p can be rescaled by using a symmetry of the 
field equations under a constant shift of accompanied by a compensating rescaling 
of the pi, and we find it convenient to normalize it as follows, 

(33) 

or equivalently set pr^ z ~ l = 1. To determine the parameters of the Lifshitz back- 
ground we insert this into the remaining field equations and set p2 = 0. From the 
dilaton equation (1521) we obtain 



2(z + 2)v^T=-^P(^J , (34) 



which requires 



Ai = ~^==, Pi = ^2(^-1)^ + 2), (35) 

and equation (13T1) then relates the cosmological constant to the dynamical critical 
exponent, 

A=-i(* + l)(* + 2). (36) 

We note that Ai is negative. The sign of Aj in (1221) determines how the coupling 
of the corresponding gauge field changes with radial location and with the dilaton 
profile in (1331) we find that F$ is strongly coupled in the asymptotic r — > oo region. 
This makes the model less satisfactory, but it is not a serious problem as long as 
Ffi) does not couple directly to anything outside the gravitational sector. Its only 
role should be to modify the asymptotic behavior of the metric from AdS to Lifshitz 
through the gravitational back reaction to its background value. For applications, 
such as holographic superconductors or fermion spectral densities, any additional 
matter fields in the model should only be charged under F$ and not under F$ . 
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3.2.2 Charged black branes in EDM model 



The charged black brane solutions in [TT] have / = g and the same dilaton profile 
fl33|) as in the purely Lifshitz background. The Einstein equation fl3Tj) then reduces 
to 

rg' + (z + 2)(g-l) = -^(^-) , (37) 

while the dilaton equation f )32|) is satisfied if 

rg+(z + 2)(g-l) = -jj=[ T ) . (38) 



It follows that A2 = V 2 — 1 an d the solution for a charged black brane with an event 
horizon at r = r is given by 

/<-> M)(r 2 +sen w 

F« = ^-1)^ + 2)^-' (^V + , (40) 
F<? = ^ r >H" + \ ° (41) 



r 

2 v / i 3 T 



(42) 



This solution is valid for any value of z > 1 and can be generalized in a straightforward 
fashion to more general spacetime dimensions [11 J. In the z — > 1 limit, the metric and 
the Maxwell gauge field F rf reduce to those of a standard AdS-Reissner- Nordstrom 
black brane, while the auxilliary gauge field F$ vanishes and the dilaton becomes 
independent of r. 

The scale transformation (fT8j) is a symmetry of the black brane metric ( 139]) . pro- 
vided ro is rescaled in the same way as r. This is a special feature of black branes 
with a planar horizon and it implies that the radial location of the horizon ro does 
not have any physical meaning by itself. The same is true for planar black branes in 
the EMP model. It is therefore important to always use scale invariant combinations 
of physical quantities when presenting results obtained in these holographic models. 



3.2.3 Asymptotic behavior 

A distinctive feature of numerical black brane solutions in the EMP model is the 
presence of a mode that is not well behaved in the r — > 00 asymptotic region [TT?|l2ip22] . 
Depending on the value of the dynamical critical exponent, it can either be a growing 
mode that takes the system away from the Lifshitz fixed point or a decaying mode 
that only converges very slowly to the asymptotic fixed point. The solutions are 
usually obtained by numerically integrating the field equations of the model, starting 
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from initial conditions near the horizon and proceeding out towards the asymptotic 
region. The initial data at the horizon is then fine tuned to remove the mode in 
question. This works reasonably well over a range of black brane temperature and 
charge density but the required fine tuning becomes progressively more difficult as 
the temperature is lowered and the black brane charge approaches its extremal value. 



EDM eigenvalues 



EMP eigenvalues 




Figure 1: A comparison of mode eigenvalues of the linearized field equations around 
a Lifshitz background in the EDM model (left) and the EMP model (right). 



In the large r limit, the analytic black brane solutions of the EDM model have 
modes that go as r~ z ~ 2 and r~ 2z ~ 2 , and the accompanying coefficients are the brane 
energy density and the charge density respectively. It turns out, however, that the 
Lifshitz fixed point is always unstable in this model and that a generic perturbation 
around it will include a growing mode that is not present in the black brane solutions. 
To see this, we linearize the field equations ( 130]) - ( 132|) around the Lifshitz fixed point. 
The function /(r) can be eliminated between the two Einstein equations, leaving a 
system involving only g(r) and <f)(r). We allow for a dilaton that varies from its 
background value, 

' r 



Mr) = { 1_ 

r 



I z— 1 f(r) 



(43) 



and write g(r) pa 1 + 7(7"). Working to first order in 7(7"), <f(r), and x(r) = r<//(r), 
we obtain 



d 


' 7 " 




dr 








. X . 





-z-2 -(z-l)(z+2) -z+1 
1 

2(z+l)(z+2) -z-2 





' 7 " 




f 




. X . 



Pi 

4 



r ) 



2z+2 



1 






(44) 



At large r, a solution to the full non-linear system will approach a linear combination 
of the eigenmodes of the linearized system plus a universal mode coming from the 
source term. The eigenvalues, 



-z-2, X - (-z - 2 ± v / (^ + 2)(9z + 10)) , 



(45) 
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are plotted on the left in Figure [TJ showing that the system has a growing mode for 
any value of z. The result of the corresponding mode analysis of the EMP model ^ZT\ 
is shown on the right in the figure, for comparison. In both models, turning on a 
growing mode would change the asymptotic behavior of the metric and take the dual 
field theory away from the Lifshitz UV fixed point. In this paper we are interested in 
models with Lifshitz scaling and we therefore only consider configurations where the 
growing mode is not turned on. 

Both models have an eigenmode that falls off as r~ z ~ 2 and from a general anal- 
ysis of holographic renormalization in asymptotically Lifshitz spacetime [23Tf27] one 
concludes that the coefficient in front of such a mode is proportional to the energy 
density of the configuration. In both models the brane electric charge density enters 
the linearized system as a universal mode due to a source term but the falloff with 
r is different. In the EMP model the charge mode falls off as r -4 , independent of 
z, while in the EDM model the falloff is r~ 2z ~ 2 and at leading order this mode only 
enters in the metric function and not in the dilaton field. 



3.2.4 Extremal Lifshitz black branes 

At the extremal value of the charge density, P2 = ±v/4(^ + 2), the metric function 
( ]39j) in the EDM model has a double zero at the horizon, 

f(r)*(z + l)(z + 2)(j--l\ (46) 

and the near horizon metric can be written in AdS2 x R2 form, 

1 



ds 
where 



(z + l)(z + 2) 



9 ,~> dv 2 



v 2 



+ dx 2 + dy 2 , (47) 



r 

v— 1, t = (z + l)(z + 2)rn t, x = r x, and y = r y. (48) 

^0 

This means that the geometry of extremal Lifshitz black branes in the EDM model is 
qualitatively similar to that of an extremal AdS-Reissner-Nordstrom black brane at 
z — 1. It follows that various holographic computations that have been performed in 
AdS-RN backgrounds will carry over to the EDM model in a straightforward fashion. 
In particular, the single fermion spectral function for probe fermions can be evaluated 
using similar techniques as have been used in asymptotically AdS spacetime [28Tf30] 
and we expect qualitatively similar results, including signals of non Fermi liquid 
behavior, to be obtained in the EDM modeo 

Less is known about extremal black branes in the EMP model. The commonly 
used method for finding numerical black brane solutions breaks down before the 



3 For calculations of fermion correlation functions in Lifshitz backrounds see [31-33 . 
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extremal limit is reached and the known exact solution at z = 4 is non-extremal. An 
analysis of the field equations of the EMP model at extremal charge density suggests 
that the near horizon geometry of an extremal black brane is in fact AdS2 X R2 in this 
model as well. The near horizon metric can be brought into the same form as fj47|) . 
but in this case the required coordinate transformation turns out to be non-analytic 
at r = r and the Proca field also goes to zero in a non-analytic fashion at r = r . 

4 Black brane thermodynamics 

The Hawking temperature of an asymptotically Lifshitz black brane is determined in 
the usual way by considering the near horizon behavior of the metric, 



,2+1 



T 



For the analytic black brane solutions of the EDM model one finds 



T 



z + 2 



(49) 



(50) 



while, as usual, only numerical results are available for the EMP model. Figure [2] 
compares the Hawking temperatures of z = 2 black branes in the two models. Other 
values of z give rise to qualitatively similar graphs. 




Figure 2: Hawking temperature as a function of the dimensionless charge density 
for z = 2 charged black branes. On the left, the dashed (blue) curve shows the 
temperature in the EDM model (equation ( 1501) ) while the solid (red) curve is from 
numerical solutions of the EMP model. On the right, the scale invariant combination 
T/p is shown for the same z = 2 data. 



The EDM action fT22|) does not depend explicitly on the gauge potential A t but 
only its radial derivative. The associated radially conserved quantity is the charge 
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density carried by the black brane, 



<5Sedm 

r" 



*-\H<**F$\ (51) 



' 5d T A? ] V / 

and inserting and from (1411) and (1421) gives 

P = P2Tq- (52) 



In the following, we have chosen to keep fixed the physical charge density p when 
calculating thermodynamic quantities and express our results in terms of scale invari- 
ant combinations such as T/p 2//2 . An alternative choice is to instead keep fixed the 
chemical potential of the Maxwell field. 




Figure 3: The ratio of specific heat to temperature in the EDM model for dynamical 
scaling exponents z = 1 (black solid), z = 2 (blue dot-dashed) and z = 3 (red dashed). 



The specific heat at fixed transverse volume and charge density is easily calculated 
in the EDM model. It can then be compared to numerical results from the EMP 
model, obtained in [31], and one finds that the two models predict different low- 
temperature behavior for this observable. The specific heat at fixed volume and 
charge density in the dual boundary theory is given by 

Ov* = T§, (53) 

where S = nr^ is the Bekenstein-Hawking entropy density of the black brane and the 
temperature T in the boundary theory is identified with the Hawking temperature 
( 1491) . Keeping p in ( 1521) fixed implies 

d d 2p 2 d 

dT = Wo ~ -WW^ (54) 
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and one then finds 

Cv, P dS/ dr 
T dT/dr 

32tt 2 / p 

4z(z + 2)-{z- A)p 2 \J 2 

Figure [3] plots the scale invariant combination p( z ~ 2 ^ 2 Cv, P /T against T j p z l 2 for three 
different values of z. In all cases the C/T ratio goes to a constant at low temperature. 
The low temperature limit can also be seen directly by going to the extremal limit 
p\ — > 4(z + 2) in the above expression for C/T, 

p ( Z -2)/2^E ->2*/ 2 (z + 2)(*- 6 >/V as T^O. (56) 

The low temperature behavior of the specific heat is very different in the EMP model. 
Numerical calculations reported in [34J show a C/T ratio at fixed charge density that 
grows as the temperature is lowered and diverges in the T — > limit for all z > 1. 

At high temperatures the specific heat in fl55l) satisfies a simple scaling law con- 
sistent with the underlying Lifshitz symmetry, 

C V , P ~ T 2 ' z . (57) 

This scaling is also seen in the EMP model [31] It follows on general grounds 
from the statistical mechanics of a system with u ~ k z dispersion in two spatial 
dimensions [35J. 




5 Two point functions at finite charge density 

We now turn our attention to finite temperature correlation functions. For earlier 
studies of scalar correlators in Lifshitz black hole backgrounds see [36] 137]. In this 
section we study the temperature dependence of two point correlation functions of 
scalar operators at fixed charge density. We first use the geodesic approximation, 
valid for large scaling dimensions, to simplify the calculation and get a sense of the 
finite temperature behavior Jf| We then consider general scaling dimensions and solve 
the scalar wave equation in the background spacetime of a charged Lifshitz black 
brane. 

4 The specific heat only depends on the geometry (through the entropy) and not directly on the 
dilaton held. More general observables, such as correlation functions of fields coupled to the dilaton, 
can fail to be Lifshitz symmetric because the dilaton itself transforms under a Lifshitz rescaling. 

5 The finite temperature two point function at vanishing charge density was obtained in the 
geodesic approximation in |20) . 
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5.1 Geodesic approximation 

First we consider operators with large scaling dimensions A. On the gravitational 
side this maps to large particle mass, for which we can approximate m ~ A. We can 
use the geodesic approximation to calculate the two point function [38J, 

(0(x)0(x')) ^ e - 2A e- A W^1i^ ; (58) 

where x M (r) is the geodesic of minimal length and e is an infrared cutoff in the bulk 
spacetime r < 1/e. 

The metrics of interest have the form fj28|) and we can parametrize the geodesies 
as r = r(x), t = t(x) and y = y(x), where x and y are transverse coordinates. 
Because of translational symmetry in the t and y directions the momenta pt oc t'(x) 
and p y oc y'(x) are conserved. By a combination of translation and rotation in the 
x-y plane we canalways arrange both endpoints of the geodesic to be at y — 0. It 
then follows that y'(x) = identically, since, if the derivative did not vanish, the 
geodesic could not start and end at the same value of y. Similarly, for an equal time 
commutator, the constant value of t'(x) must in fact be zero for otherwise the two 
endpoints would be at different values of t. 

In this case, the geodesic length functional simplifies, 



S = A / dxJr(x) 2 + -¥^¥- = A I dxL. (59) 
\i r{x) 2 g(r) J 

Since the "Lagrangian" L has no explicit x dependence, there is a conserved Hamil- 
tonian, 

H = - L = -4- (60) 
dr' L K J 

The value of the Hamiltonian can be calculated at the turning point of the geodesic 

at r = r*. By translation symmetry the turning point can be taken to be at x = 

and then we can use r'(0) = to give H = — r*. This way we are lead to a first order 

differential equation for the geodesic, 

( r ')2 = ^ f^-lY (61) 



The on shell particle action is now given by 



L 2 A 

S = A / dr- = — / dr = (62) 



where we have used the chain rule and the equation of motion (!6T|) . It is straight- 
forward to solve (ED) numerically and evaluate the integral (1621) . Before doing this, 
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we note that the two point function can be computed analytically for large distances. 
First we write (loTj) in integral form 

i , r' 1 1 

— |x — x'| = / dr — -j^^^^^^^. (63) 



2 



gr 



At r 



2 



The key is to note that large distance |x — x'| corresponds to a geodesic with a turning 
point approaching the horizon r = r§. In the limit r* — > r$ the integrals ( 1621) and (1631) 



both diverge logarithmically at the lower endE So in this limit we can approximate 

i, ,. r" 1 , i i r , i 

x — x|= / ar — j pa — / ar — ; = + (64) 



2 "" < 9 r^-i) •- ^{j-i, 

where the dots denote terms that are subleading as r* — > To- This way we get 

S = Ar |x-x'| + ... (65) 

where the dots denote terms that increase slower than linearly in the limit of large 
|x — x'|. Using (158j) we obtain 

G(x, x') = (£>(x, t)0(x!, t)) oc e-l x - x 'l/«, (66) 

where the correlation length is given by 

5 = s;- (67) 

In order to obtain the temperature dependence of the correlation length we need to 
know the relation between the temperature and the position of the horizon. This 
has to be obtained mostly numerically. The temperature is related to the position 
of the horizon through (H9j) . keeping p in (152|) fixed. Figure H] shows the results of a 
numerical evaluation of (|67|) for both holographic models function of T/p 2 / 2 . A 
numerical evaluation of the full correlation function in the geodesic approximation for 
z = 2 is shown in Figure Correlation functions for other values of z and different 
values of T j p z l 2 are qualitatively very similar to the ones shown in the figure. As can 
be seen from Figure the numerical solution indeed agrees with our (approximate) 
analytic result fIBTj) at large distances. 

The short distance behavior of the correlation function has the form 



6 The integral in (|62|) is also divergent at the upper limit when the cutoff is sent to zero. This 
divergence is independent of r* and has the same value in the vacuum. We can ignore it as it cancels 



with the explicit power of e in ([58 
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Figure 4: Correlation lengths in the two models as a function of T/p 2 ^ z for the 
dynamical critical exponent z = 2 (top) and z = 3 (bottom). The solid red curves 
correspond to the EMP model and the dashed blue curve corresponds to the EDM 
model. 



independently of the temperature. 

As can be seen from Figure HJ the correlation length approaches a finite non- 
vanishing constant value as T — > while keeping the charge density fixed. For a 
vanishing chemical potential the correlation length diverges as £ oc T- l / z as Figure H] 
indicates. 

5.2 General scaling dimensions 

When the scaling dimensions of the operators in the two point function are small, the 
geodesic approximation no longer applies and one has to solve the full wave-equation. 
The Euclidean action for the boson is taken to have the standard form@ 



c/ 3 xv/^70 2 | r=e -i, (69) 



7 We note that the counter term in (|69|) is not the most general one, but it applies to the values 
of the scaling dimension A < 2 + z/2, to which we specialize in what follows. More generally one 
has terms including derivatives with respect to the transverse dimensions [5]. 
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A _l Log(G) A 'log(Go) 




Figure 5: Logarithm of the two point equal time correlation function as a function 
of distance for z = 2. The dashed line corresponds to a reference line with slope -1. 
The left figure is in the EDM model and T j p z l 2 rs 0.0229 while the right figure is in 
the EMP model with T/p z / 2 « 0.109. Correlators for other temperature values are 
almost indistinguishable from the ones shown in the figure. Most of the dependence 
on z and T comes through £. 



where 7 is the determinant of the induced metric at r = e 1 and A_ is given below 
equation fl75l) . The two point function is obtained by solving the wave equation, 

1 



-9 



-MV=gg ,tt 'd v <i>)-m i <j> = o, 



(70) 



in the bulk spacetime. We are interested in correlation functions at thermal equilib- 
rium, and thus we work in Euclidean signature. This means in particular that the 
field 6 should be Fourier transformed as 



E 



d 2 k 

(2tt)' 



-iu; n T+ik-x 



</>n(r, k), 



(71) 



where u n = 2irTn are the Matsubara mode frequencies. To obtain the bulk to 
boundary propagator we require regularity at the horizon. For the numerics we find 
it convenient to use the coordinate u = 1/r in terms of which the metric reads 

-2 1 , J„.2 



, 2 „. s dr 2 If du 2 , 2 \ 
ds 2 = f( u )— + — — - + dx 2 . 

u lz u l \q[u) / 



u 2 \g(u) 

We are lead to solve the wave equation 

J 



u 



3+z 



— 2z 

jdviu-^y/fjdvfa) - (^col + u 2 k 2 + m 2 )0 n = 0, 



(72) 



(73) 



together with the regularity condition, that, close to the horizon u = u = l/r , the 
field behaves as 



(f) n (u, k) w exp 



-ul 1 uj n 



du' 



y/f{u')g{u') 



(74) 
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Figure 6: The momentum space correlator for z = 2 and u n = in the EMD model 
at finite temperature and zero charge density. The scaling dimension of the scalar 
operator is A = 2 + y / 3/2. 



where has a regular Taylor expansion in u — Uq. The modes 4> n (u, k) have the 
standard asymptotic behavior as u — > 



(75) 



where A± = (2 + z)/2 ± v and v = a/(2 + z) 2 /4 + m 2 . The holographic dictionary 
with our Fourier transform conventions reads 



e " 5e = ( exp 



n ^ 

which leads to the two point function 

(O n {k)O n ,(k')) = 



5 2 S F 



where 



p 2 s^-m^-k') 

T8 n+n ,, 8 2 (k + k')(2n) 2 G n (k), 



(76) 



(77) 



G n (k) = 2v^- y 



(7S 



For the alternative quantization P9"|HU] (which corresponds to treating the subleading 
mode (f)^ as the source) we have 



G n (k) = -(2v) 



6 H 

-1 'rn 



(79) 
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The momentum space correlator has poles on the imaginary k axis, which means that 
the real space two point function decays exponentially, G(x, x') oc e - fe *l x - x 'l ) where 
A;* is the pole closest to the real k axis. It is sufficient to take into account the lowest 
Matsubara mode as the higher ones have a faster fall of in position spacfl We solve 
( 175]) numerically and extract the momentum space two point correlator using fl78j) . 
with n = 0. Results from a numerical evaluation of the finite temperature two point 
correlator is shown in Figure El This indeed shows poles on the imaginary k axis. 
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Figure 7: The correlation length, as extracted from the position of the lowest pole of 
the momentum space correlator, as a function of T/ p z l 2 . The dashed (blue) curves 
are in the EMD model while the solid (red) curves are in the EMP model. All the 
plots are for z = 2 and A = 2 + y/3/2. 



By tracking the position of the lowest pole in the momentum space propagator, 
we can obtain the correlation length as a function of T/p z l 2 as the temperature is 
varied. The correlation length for z = 2 obtained this way is shown in Figure [7J for 
both holographic models. 

6 Double trace deformations 

In the previous section we considered the finite temperature dependence of scalar 
correlation functions at finite charge density. While working at finite charge density 
gives rise to non-trivial dependence on the dimensionless ratio of the length scales 
associated to p and T, it has the drawback that at zero temperature there is still a 
finite length scale in the problem, associated with the charge density. 

Another way to study the temperature dependence of the correlation length is 
to work at zero charge density and add a double trace deformation for the scalar 

8 The same procedure has been used earlier in calculating the Debye screening length using holog- 
raphy in JV = 2* supersymmetric Yang-Mills theory |41) . 
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operator |421I43]. 

§S E = I I d 6 xXO{xf. (80) 



2 

The coupling constant A introduces a new reference length scale into the problem, 
which, when combined with T, gives rise to a dimensionless ratio that enters into ther- 
mal correlators. For a concise review of double trace deformations see the Appendix 

of m. 

A prescription for calculating two point correlation functions in AdS / CFT in the 
presence of a double trace deformation was given in [45J. The same argument goes 
through in asymptotically Lifshitz spacetime and we do not repeat it here. Denoting 
the original momentum space two point function as G n (k), the two point function in 
the deformed theory is found to be 

G (X) {k y G n(k) 



1 + XG n (k)' 



The form of the two point function ( IHTj) can also be understood in a simple way from 
large- iV factorization of correlation functions. The deformed two point function is 

O n (k)O n ik')e-^ J ^ XOn{q) °-" { - q) ) = T6 n+n ,, 6 2 (k + k')(2n) 2 G£\k), (82) 



which after expanding the exponential in a power series, and using large-N factoriza- 
tion, becomes 

771=1 V ' 

recovering (IHTj) . 

The method that was described in the previous section for solving for the two point 
function G(k) applies here as well. Since we are interested in relevant deformations of 
the theory, we consider correlation functions in the alternative quantization [3^||4D]. 
for which the undeformed correlator is given by f!79l) . In this case the deformed 
correlator becomes 

^ {+) (£;)\-i 



G ^-^m> • <84) 



By tracking the lowest pole of the deformed two point function (J84J) we obtain the 
correlation length shown in Figure |8j At high temperature the dimensionless corre- 
lation length ^T 1 / 2 is simply that of the alternative quantization without a double 
trace deformation. As the temperature is lowered, the dimensionless combination 
^T 1 / 2 flows from its value in the alternative quantization, to its value in the stan- 
dard quantization. This reflects the renormalization group flow from the alternative 
quantization to the standard one [4D]|44]. In the vacuum the renormalization group 
flow is seen when considering the behavior of the correlation function as a function 
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Figure 8: The correlation length as extracted from the position of the lowest pole of 
the momentum space correlator in the double trace deformed theory. The solid red 
curves are in the EMP theory while the dashed blue curves are in the EDM theory. 
Both plots are for z = 2 and A = 3/2. 

of distance, short distances corresponding to the alternative quantization and large 
distances to the standard quantization. In our case the renormalization group flow is 
instead revealed SIS db function of the temperature. 

7 Discussion 

In this paper we have considered two bottom up holographic models both of which 
admit a metric with a Lifshitz scaling symmetry, with arbitrary dynamical critical 
exponent, and include a U(l) gauge field. Both models have charged black brane 
solutions that give access to finite temperature physics in the corresponding dual field 
theories. At high temperature (compared to the scale set the U(l) charge density) the 
behavior of thermodynamic quantities is governed by the underlying Lifshitz scaling 
symmetry. As the temperature is lowered, keeping the charge density fixed, various 
thermodynamic quantities in the two models start to differ from each other. In 
particular, a sharp distinction between the two models appears in the specific heat, 
which for the EMD model behaves at low temperatures as 

C 

— — > const, as T — > 0, (85) 

while numerical results suggest that C/T grows without bound in the limit of zero 
temperature in the EMP model [51] . It thus appears that the same relevant deforma- 
tion, by a chemical potential, takes us to a different kind of a fixed point in the zero 
temperature limit in the two models, at least as far as thermodynamics is concerned. 
The infrared (near horizon) geometry of extremal black branes has an AdS2 x R 2 
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region in both of the models, as discussed in section [3j but in the EMP model only 
after a singular coordinate transformation. 

For the finite temperature equal time two point function, at finite charge density, 
we find the standard power law behavior at short distances, while for large distances 
it decays exponentially due to the presence of a horizon in the bulk. The corre- 
lation length extracted from the exponential decay is shown in Figure H] for large 
scaling dimension operators and in Figure [7] for generic scaling dimensions. For high 
temperatures T j p z l 2 > 0.2, the temperature dependence of the correlation length is 
consistent with scale invariance, implying a large temperature range where the critical 
point dominates the physics. As expected, the two models agree quite well in this 
range. On the other hand, the correlation length exhibits distinctly different behavior 
in the two models at low temperature and fixed charge density. In the EMD model 
we obtain the following behavior in the geodesic approximation, 



while the same quantity diverges in the zero temperature limit in the EMP model, 



In both cases £ tends to a finite value as T — > when the charge density is kept fixed 
at a non-zero value. 

The correlation length in the presence of a double trace deformation is shown 
in Figure [H In this case the correlation length diverges in both models in the zero 
temperature limit, while keeping the double trace coupling A fixed. The temperature 
dependence of the correlation length has a scale invariant form in two regions, at high 
temperatures and at low temperatures. The combination ^T 1 ^ runs between its fixed 
point values, from the alternative quantization at high temperatures to the standard 
quantization at low temperatures. Overall, the double trace deformation leads to a 
qualitatively very similar behavior in both holographic models. 

The double trace deformation leads to a Lifshitz symmetric result in the zero tem- 
perature limit while at fixed charge density there remains a finite length scale in the 
problem even at zero temperature. We can therefore more meaningfully compare the 
double trace deformation case to the quantum Lifshitz model, which, as was reviewed 
in section 121 has several temperature dependent length scales. The scale associated 
with a vortex plasma in the quantum Lifshitz model behaves as t; VO rtex oc e~ Ec ^ 2T , 
which changes more rapidly with T than we find in the holographic correlation func- 
tions. The comparison to the "intermediate" length scale £r in (I14p. given by 




(86) 
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works somewhat better. The dimensionless combination ^T 1 / 2 is seen to be a fairly 
mild function of temperature, even if it diverges in the zero temperature limit. At 
a qualitative level, the double trace deformation gives rise to similar behavior as it 
induces a mild temperature dependence on ^T 1 / 2 . 

Both types of deformations in both holographic models show similar behavior for 
the combination ^T 1 / 2 . It is seen to be an increasing function of temperature in all of 
the cases we have studied here. It would be interesting to understand this from the 
perspective of the finite temperature renormalization group. 
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